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Abstract 

By using relations derived from renormalization group equations (RGEs), we 
find that strong indirect constraints can be placed on the top squark mixing 
phase in At from the electric dipole moment of the neutron (d n ). Since mt is 
large, any GUT-scale phase in At feeds into other weak scale phases through 
RGEs, which in turn contribute to d n . Thus CP-violating effects due to 
a weak-scale At are strongly constrained. We find that |Iim4^ w | must be 
smaller than or of order llmP^^I, making the electric dipole moment of the 
top quark unobservably small in most models. Quantitative estimates of the 
contributions to d n from A u , A^ and B show that substantial fine-tuning is 
still required to satisfy the experimental bound on d n . While the low energy 
phases of the A's are not as strongly constrained as the phase of B EW , we 
note that the phase of a universal A GUT induces large contributions in the 
phase of B EW through RGEs, and is thus still strongly constrained in most 
models with squark masses below a TeV. 
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I. INTRODUCTION 



Supersymmetry (SUSY) |T| is one of the most compelling extensions of the Standard 
Model. It is the only known perturbative solution to the naturalness problem B], it unifies 
the gauge coupling constants for the observed value of sin 2 6^, it allows radiative EW 
symmetry breaking, and the lightest SUSY partner provides a good dark matter candidate. 
SUSY models with such features are generally in excellent agreement with experiment, and 
there is even the possibility that a recent CDF event || is of supersymmetric origin Bj]. 

One of the few phenomenological problems associated with SUSY models is their generi- 
cally large predictions for the electric dipole moment (EDM) of the neutron, d n . Supersym- 
metric models with universal soft breaking parameters have two physical phases, beyond 
the CKM and strong phases of the SM, which can be taken to be the triscalar and bis- 
calar soft breaking parameters A and B. These phases give a large contribution to d n , of 
order lCr 22 (100 GeV/M susy ) 2 e cm, where M susy is a characteristic superpartner mass. The 
experimental upper bound on d n is of order 10 _25 e cm ||, so that if superpartner masses 
are near the weak scale, the phases of these complex soft parameters must be fine-tuned to 
be less than or of order 10~ 2 -1CU 3 since there is no a priori reason for them to be small 
0. If one wants to avoid such a fine-tuning, there are two approaches: suppress d n with 
very large squark masses (greater than a TeV) [0], or construct models in which the new 
SUSY phases naturally vanish ||. Models with very heavy squarks are unappealing because 
in such models LSP annihilation is usually suppressed enough so that the relic density is 
unacceptably large ||. They also lead to a fine-tuning problem of their own in getting the 
Z boson mass to come out right in EW symmetry breaking. 

It is natural to consider solutions of the second type, and demand that the soft phases 
are zero by some symmetry. While that would leave only a small CKM contribution to 
d n pH|-p^|, and thus avoid any fine-tuning in meeting the experimental bound on d n , it 
would also mean that there is no non-SM CP violation, which is needed by most schemes 
for electroweak baryogenesis fll4 |. Also, such models do not generate signals of non-SM 
CP violation, such as those involving top squark mixing. There are ways of naturally 
obtaining small nonzero soft phases which leave sufficient CP violation for baryogenesis 
|T5Hl8|, but these phases would still have to meet the bounds from d n and would probably 
be unobservably small in most EW processes — unless the soft terms are not universal. 

Recently it has been pointed out that large non-SM CP- violating top quark couplings 
could be probed at high energy colliders fL9 |. A measurement of a large top quark EDM, for 
example, would indicate physics beyond the SM, and it is interesting to ask whether SUSY 
models can yield an observable effect. Several references have attempted to use CP violation 
from top squark mixing due to the complex parameter A t to yield large CP- violating effects 
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in collider processes involving top quarks p0 |. Such papers either explicitly or implicitly 
assume nonuniversal soft couplings A q at the GUT scale; otherwise, the phase of A t would 
be trivially constrained by d n . We consider whether it is possible to obtain large effects due 
to the phase of A t at the EW scale by relaxing the universality of A. We will show that 
due to renormalization group induced effects on other low energy phases, the phase of A t 
is strongly constrained by d n , and it is not possible, for most areas of parameter space, to 
have large CP- violating effects due to the imaginary part of A t . 

We will assume that no parameters are fine-tuned and thus we will require the phases at 
the GUT scale to be either identically zero (presumably through some symmetry) or no less 
than 1/10. If one permits an arbitrary degree of fine-tuning, the whole SUSY CP violation 
issue becomes moot, and one can derive no constraints on the phase of A t . While one can 
construct models which give small universal phases, as we said above, the fine-tuning needed 
to evade the constraints we derive is unlikely to be explained naturally. Our approach in 
this paper is to assume the reasonable fine-tuning criterion we have just outlined, and ask 
what it implies about low energy SUSY CP- violating phenomenology. 

In Sec. II, we review the basics of SUSY CP violation. We present our results derived 
from RGEs in Sec. Ill, and impose the neutron EDM constraints on lmA t using those results 
in Sec. IV. In Sec. V we discuss top squark mixing induced CP violating observables in 
more detail in light of our constraints on the phase of A t , and we give some concluding 
remarks in Sec. VI. The details from Sec. Ill are written up in Appendix A, and the full 
one-loop calculation for the SUSY contribution to the neutron EDM is given in Appendix 
B. 



II. SUSY CP- VIOLATING PHASES 

The soft breaking potential in the MSSM is 
- C so ft = - \mi\ 2 \ipi\ 2 + - M A AA + 

A 

e^AuU^YuQ^Hl + e l3 \A D D\Y^ L \lP d + e l3 \A E E R YlV L \H\ + (1) 
tijB^iHlHi + h.c. 

(2) 

where we take A v = diag{A u , A c , A t }, A D = diag{A d , A s , A b }, A E = diag{A e , A^, A T }\ Y n , 
Yd, and Ye are the Yukawa coupling matrices; Q, L, Ur, Dr and Er are the squark and 
slepton fields; A are the gauginos and <pi are the scalars in the theory. 

A common simplifying assumption is that this soft Lagrangian arises as the result of a 
GUT-scale supergravity (SUGRA) model with universal soft triscalar coupling A, gaugino 
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mass M\ = M\m, and scalar mass rrii = mo. This provides an explanation for the absence 
of flavor changing neutral currents which arise from loops with squarks of nondegenerate 
mass ||21|| . Such supersymmetric models have only two independent physical CP- violating 
phases beyond the CKM and strong phases of the SM |T(J although these phases appear in 
several different linear combinations in low energy phenomenology []I7||22]. We will take the 
two physical phases to be ArgA and ArgB. 

It turns out that all CP violating vertices in this model arise through the diagonalization 



of complex mass matrices [|15|]. The complex quantities which appear in these matrices are 
A q + fi*R q and //, where R q is tan/3 (the ratio of Higgs vacuum expectation values) for 
q = d, s, b and cot (3 for q = u, c, t, and where the phase of /i is simply equal to the 
phase of B* by a redefinition of fields. Thus for d n , which involves only u and d quarks, 
there are only contributions from three low energy combinations of the two SUSY GUT 
phases: Arg(A d — /itan/3), Arg(A M — //cot/3), and Arg/i. (In the Appendix B, a complete 
expression of d n is given which includes suppressed contributions from phases of the other 
squark mixings.) 

Even with universal boundary conditions, the elements of the matrices Ajj, Ad and 
Ae have distinct phases at the EW scale because of renormalization group evolution. We 
will also relax, in some places, the assumption that their phases started the same at the 
GUT scale. We assume (for simplicity) that these matrices are diagonal. One possible 
consequence of this approach is that one could have d n ~ because ImA^ and Iim4 n ~ 0, but 
other A q , notably A t , could have large phases which lead to observable effects. These include 
angular correlations and polarizations [pO] , including effects attributable to the electric dipole 
moment of the top quark, d t . As discussed in the Introduction, this scenario is strongly 
constrained by RGE running. 



III. RENORMALIZATION GROUP FLOW OF COMPLEX SOFT TERMS 

The goal of this section is to demonstrate how a large phase in A t can feed into other 
parameters in the theory through renormalization group running. The imaginary part of A t 
at the weak scale, Af w , is determined by running Af UT (and for large tan/3, Af UT ) down 
to the weak scale via the renormalization group equations (RGEs). (For compactness of 
notation, we will define x = Imx in the following sections.) We will show that large Af w 
induces potentially large values of B EW and A^Y, which give an unacceptably large neutron 
electric dipole moment. 

Rather than write RGEs for the whole effective theory, we need only consider a complete 
subset of them which includes A q and B. The running of these soft terms depends upon 



3 



the gaugino masses, the top and bottom Yukawas (we ignore tiny effects from the other 
Yukawa couplings) and the gauge coupling constants a a = A^/47r (a = 1, 2, 3). We define 
t = l/47rln(Q/M GC/T ) and write 

dM a 



dt 
dA 

dt 

dA u ^ c 

dt 
dA b 

dt 

dA djS 

dt 
dB_ 

~dt 

dott 

~dt 

dab 

~dT 

da a 



2b a a a M a (3) 

2c a a a M a + Yla t A t + 2a b A b (4) 

2c a a a M a + 6a t A t (5) 

2c' a a a M a + 2a t A t + 12a b A b (6) 

2c' a a a M a + 6a b A b (7) 

2c> a M a + 6a t A t + 6a b A b (8) 

2a t {-c a a a + 6a t + a b ) (9) 

2a b (-c' a a a + a t + 6a b ) (10) 



dt 

where a is summed from 1 to 3, and 



2b a a 2 a (11) 



&«=(— ,1,-3), (12) 

(13) 



H i 

5 ' ' 


-3 
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ci=^,3, T j, (14) 
<'=(l3,0), (15) 



and the Yukawa coupling constants a tjb = Af 6 /47T are related to the masses by 

x 92 m 1 5-2 m fe 1 

A* = —7= 77, A& = — ^= -. (16) 

V 2 mvi/ sin p \f 2 m w cos p 

We note that some references J23| list the atA t coefficient in ([5]) as 2, but we have confi- 
dence that the coefficient is actually 6 [13,23]. Nevertheless our conclusions do not depend 
qualitatively on this coefficient. 

We are mainly interested in the evolution of A q and B. We can set the phase of the 
common gaugino mass to zero at the GUT scale by a phase rotation and then Mj = at all 
scales. Therefore the RGE for the imaginary parts of the A q and B can be written without 
the M a terms: 
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^ = 12a t A t + 2a b A b , (17) 
dA b 



2a t A t + 12a b A b , (18) 
6a t A t , (19) 
Qa b A b , (20) 



dt 
dA u 

~dT 

dA d 

dt_ 
dB 

— = Qa t A t + Qa b A b . (21) 

Using the above RGEs, we can derive the following useful relations: 

AB =AA ViC + AA d)a = ^(AA t + AA b ), (22) 

AA UjC = A ( 6 AA - AA b ) , (23) 

AA d , s = A ( 6 AA - Alt) , (24) 

where AS = £? Gt/T - B EW } etc. For small tan /3, we can neglect m& so that these relations 
simplify to 

AB =AA UiC =3AA h = -AA t , 
AA dtS = 0. (25) 

Thus, given the GUT values, to obtain the low energy values for the imaginary parts of all 
the soft terms, one only needs to find Af w and Af , and for small tan /3, we only need the 
former. 

In the small tan/? limit (a b ~ 0), we can use Eq. ( |i~7|) to obtain the ratio of EW to GUT 
scale values of the imaginary part of A t : 



r t = Af w /Af UT = exp 



12a t (t)dt 

tEW 



(26) 



If the top quark were light, the integral in Eq. ( p6|) would be small and r t would be close 
to one, but since the top quark is heavy, we find that r t is well below one. We can use the 
relations in Eq. ([25]) and the definition for r t in ( p6[ ) to relate the low energy values for the 
imaginary parts of A t to B and A u (for small tan (5) : 

Af w = Z^L (b ew - B GUT ) = T ^ (A EW - A° UT ) . (27) 

We will make the simplifying assumption that A^ UT and B GUT are zero. As we will see in 
the next section, this is reasonably well justified by our fine-tuning criterion, at least for the 
phase of B. 
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Next, we must find r t . We obtain a pseudo- analytic solution to Eq. ( p6|) in terms of 
EW and GUT scale quantities in Eq. ( |39| ) of Appendix A, but this is useful only if one has 
already obtained the GUT values for the a's by numerical integration of the RGEs. While 
we cannot find a truly analytic solution to Eq. (p6|), we can place an analytic upper bound 
on r t which is sufficient to make our point. We note that the integral in Eq. ( p^) is simply 
the area under the curve of the top Yukawa a t as it runs from the EW scale to the GUT 
scale. Thus we can place an upper bound on r t simply by finding a lower bound to that 
area. In Appendix A, we do this by placing a lower bound on at(t) at each t, and we obtain 

r t < 1 - 12a EW /f EW , (28) 

which is valid for small tan/3 so long as 12a EW < /ew- Here Jew = 2c a a EW ~ 1.5 + 
32/3(af w/ — .12), so, for example, Eq. fl28|) is valid for m t = 175 if 1.3 < tan/3 <C m t /m b 
(for smaller tan/?, r t gets closer to zero, but does not actually reach it). Thus we have 
placed an analytic bound on the running of A t completely in terms of EW quantities. For 
a s (Mz) = .12, sin/? — > 1 (moderate tan/3) and m t = 175 (m t = 160), we find that r t < .43 
(r t < .52), which, from Eq. (0), corresponds to \Af w \ < l.h\B EW \ (\A EW \ < 2.2\B EW \). 
For small tan/3, the bound is even stronger, so that for tan/3 small enough to neglect nib 
effects, we obtain 

\Af w \ <2.2min {\B EW \,\A EW \} (29) 

and in practice the coefficient is less than 2. 

In Fig. 1, we plot r t (= Af w /Af UT ) as a function of the top Yukawa coupling for 
different values of a s (Mz) in the limit where effects proportional to mt, can be ignored. For 
m t > 160 GeV, X t is always greater than about 0.87 for all values of tan/3, which means 
that r t is always less than .45, in agreement with our analytic bounds. Also plotted are 
-B EW /Af UT = (1 - r t )/2, and -B EW /A EW = (1 - r t )/2r t , which is greater than 1 (0.6) 
for m t = 175 (160). Thus |v4f H/ | ^$ l-B'^l, in agreement with our analytic results. 

Next we consider moderate tan/3, where one must take into account the mixing of A t 
and At, but where tan/3 is not of order m t /mb. For Af UT /Af UT > 0, these effects lower r t , 
and one can simply use the nib = upper bound on r t derived above.[] 

For Af UT /Af UT < (recall that with universal A this ratio would simply be +1), one 
simply maximizes the positive contribution to r t from A b to obtain (see Appendix A) 

*There is a subtlety for the case of small positive Af UT /Af UT for which there can be a net 
positive contribution to rt if Ab runs down below zero. However, the maximum effect on the bound 
is very small. 
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r t < 1 - l2ar/f EW - - 6 {A^/A^) °» ■ (30) 

Note that the last term raises the upper bound on r t , but the effect is small until tan/3 gets 
quite close to m t /m b . For m t = 175 GeV, Af UT = —Af UT , and tan/3 = .7m t /m b ~ 35 
(recall that we are evaluating all quantities at the EW scale, so m b is somewhat lower than 
the value at q 2 = ml), we find the bound r t < 0.6. 

Effects due to m b are evident in Figs. 2-4, which show Af w , B EW , A EW and A EW , 
normalized to Af UT , as a function of tan /3 for various GUT-scale boundary conditions. In 
Fig. 2, only the phase of Af UT is non-zero, while in Figs. 3 and 4, A b UT has values of +Af UT 
and — Af UT respectively. In all cases, r t (the solid curve) remains below 0.35 and has its 
largest value just below tan/3 = m t /m b for Af UT /Af UT < (Fig. 4), in agreement with our 
analytic results. This means that the EW value for the phase of A t is constrained to be less 
than about a third, independent of constraints from low energy CP violating observables. 
The magnitude of the imaginary part induced into B EW / Af UT by A t is greater than 0.35 
except for large tan/3 and A b UT /Af UT < 0. For tan/3 = m t /m b and A b UT / Af UT < 0, B EW 
actually goes through zero, because AB gets equal and opposite contributions from AA t 
and AA b there. At that point the "t" and "b" RGE coefficients are almost exactly the same 
at each t (because the Yukawa coupling runnings differ only in a small U(l) coefficient), and 
the boundary conditions have opposite signs, so that A t (t) ~ — A b (t) for all t. Of course 
A EW and A EW are non-zero because they involve different linear combinations of AA t and 
AA b , so there is still a strong constraint on A t from d n there. 

Finally we note that for large tan/3, one can place constraints on A b UT as well, since 
it can then affect other low energy phases through renormalization group running. For 
tan/3 ~ m t /m b , the constraints are of the same order as on Af UT , while for small tan/3, 
Af UT is unconstrained (though the A t contribution to A EW for small tan/3 is constrained 
to be small and A EW can be large only if \Af UT \ > \Af UT \). 



IV. BOUNDS FROM THE NEUTRON EDM 

Now that we have placed an upper bound on the magnitude of A t in terms of A u , Ad, 
and B, we need to explore the constraints on the latter three imaginary parts (in low energy 
observables, we will drop the label EW). As we mentioned in the Introduction, one of the 
strongest constraints on CP violating phases is the electric dipole moment (EDM) of the 
neutron, d n . In Appendix B, we write expressions for the full supersymmetric contribution 
to d n . One sees that all the pieces are proportional to A u , Ad, or Jx (except for the negligibly 
small pieces proportional to A q ). We can redefine the Higgs fields so that the phase of /i is 
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just the opposite of the phase of B, and thus 



B 



' n A t -B GUT ), (31) 



2r t 

where the RHS follows for small tan f3. 

In order to estimate the size of Jx we will need an estimate of \fi/B\ in Eq. fl3~T|). We 
can find this ratio by considering the two equations which \x and B need to satisfy to ensure 
that EW symmetry breaking occurs and that the Z boson gets the right mass: 

2Bfi = -{m 2 Hu + m 2 Hd + 2/i 2 ) sin 2/5, (32) 



/i 2 = ^ + — ^ -. (33) 



m| m 2 ^ - m 2 ^ ta n 2 /? 
2 tan 2 /5 - 1 

In the limit that tan/3 — > oo, we see that the right hand side of Eq. ([32]) goes to zero, so 
5^0, whereas /i 2 is not forced to zero. For tan/5 — > 1, the right hand side of Eq. ([33]) blows 
up forcing /i to take on very large values. When /i 2 dominates Eq. (^) and tan (3=1 then 
we are led to a value of \B\ = \fi\. So in both the tan/5 — > oo limit and the tan/5 — > 1 limit 
we find that > We have run thousands of models numerically p5| which include 
the one-loop corrections to Eqs. ( |32"D and (p3|) and found that > |5| is indeed a good 
relationship for most of the parameter space. As expected, it is violated most strongly for 
intermediate values of tan/5. For example, for tan/5 = 10 we have found a small region of 
parameter space where is as low as 0.4, although most solutions prefer > 1. 

We will assume that > 1, and thus the fine-tuning constraint on the phase of B is 

even stronger than on what we obtain below for the phase of /i. 

From Appendix B, we see that d n can be written in terms of the three imaginary parts, 



-, (34) 
m 



10 zo ecm tuq m m m mo 

where we have normalized the RHS by the SUSY mass scale mo, and the LHS by the region 
of the experimental bound so that the coefficients k are dimensionless. We can rewrite the 
EW imaginary parts in Eq. (Rl) using Eq. (B3) as 



10 25 ecm 10 25 ecm 2rv 



— , (35) 
mo 



where d GUT /10~ 25 e cm is just Eq. ( |34]) with EW values of A u ^ and B replaced by GUT 
quantities. It vanishes if A GI d TT and B GUT are zero. In supergravity models, \A GUT \ and 
\B T \ are of order m , so that barring fine-tuned cancellations, the GUT scale phases must 
be less than order l/k n . If the k's are greater than order 10, then our fine-tuning criterion 
dictates that we set the GUT phases to zero (presumably protected by some symmetry). 
Thus we need an estimate of the k' n s. 



8 



In Figs. 5a, 5b, and 5c, we plot the values for k^ u , k^ d , and k£ respectively in many 
different models as a function of squark mass, and as a function of tan f3 in Figs. 5d, 5e, and 
5f. We see that k^ u and k^ d are fairly flat functions of tan/5, whereas — k^ increases with 
tan/5 due to the //tan/? terms in the expression for d n . We also see that most models give 
k^ > 2 (0.8), > 7 (3), and |A#| > 100(40), for squark masses below 500 GeV (1 TeV), so 
that order one phases in all the SUSY complex quantities usually give a neutron EDM which 
is of order 100 (40) times the experimental bound. We note that these are substantially larger 
contributions (and thus stronger constraints) than claimed by the recent work of Falk and 
Olive |26|] , though this is probably due to the fact that they use very heavy squark masses in 
an effort to find the smallest fine-tuning of phases consistent with cosmology. While one can 
argue whether or not the bounds on the phases of A u ^ represent a fine-tuning, the bound on 
the phase of /x (and thus B EW , which comes from B GUT and Af ) certainly does. Thus, 
by our fine-tuning criterion, the phases of B GUT and Af UT should be zero. We note that in 
the case of universal A it is irrelevant whether or not the low energy phases of A u and A^ 
are strongly constrained, since the phase of the universal A GUT makes a large contribution 
to the low energy value of /2 (since A GUT = A GUT ). 

To give an idea of what level of neutron EDM one expects with different initial as- 
sumptions, we plot in Fig. 6 d n /l0~ 25 e cm with universal |A G ' C7T | for three cases: (a,d) 
ArgAf UT = ArgA GUT = 0.1 and all other phases zero, (b,e) ArgAf UT = -ArgA GUT = 0.1 
and all other phases zero, and (c,f) universal phases ArgA GUT = AigB GUT = 0.1. As one 
can see, even with phases of order 0.1, most models have an absolute value for d n /10~ 25 e cm 
greater than one, inconsistent with the experimental bounds. 

As can be gathered by the spread of points in the scatter plots and the number of 
parameters involved, the results depend on one's model assumptions. For example, if one 
requires tan (3 to be small (say because of b-r unification), and the squarks are allowed to be 
very heavy, then there is very little fine-tuning needed for the current experimental bound 
on d n . On the other hand, if SUSY is detected at LEP 2 or TeV 33, then even the smallest 
tan f3 models would require fine-tuning. 

In minimal supergravity models the natural scale for the A terms is m$. In Fig. 7 we 
have plotted o? n /10 _25 e-cm versus ImAf w /mo to succinctly demonstrate how quickly the 
EDM rises when ImAf w ^ 0. To construct this plot we chose a random phase for A t at 
the GUT scale, forced all other phases equal to zero at that scale, and then ran all the 
parameters down to the weak scale. A sharp drop in d n occurs at ImAf w /m$ ~ because 
can be small there and thus induces only small phases into the other low energy 
soft parameters. Models with d n around 10 e-cm at ImAf w /m ~ occur for low tan (3 
where lmA GUT 3> lmAf w but where d n otherwise tends to be smaller. This means that 
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most models with d n below the experimental bound in Fig. 6 also have a small EW value for 
Inn4 t , and thus from Fig. 7 we can place a stronger constraint on \m.Af than we obtained 
on lmA? UT : ImA? w /m < 1/20. 

Thus we conclude that models with universal GUT-scale phases of the soft parameters, 
and models in which only Af UT has a non-zero phase, have difficulty meeting the bounds 
from d n and our fine-tuning criterion. Models with non-zero 

A GUT^ A GUT^ or A GUT 

can 

meet the constraint from d n without fine-tuning, as can those with non-zero A^ UT for small 
tan f3. For the remainder of the paper, we will for simplicity set all the GUT-scale phases to 
zero except for that of A t . Even though our fine-tuning criterion implies that Af UT should 
be zero, we find it useful to ask what effects one would have if one allows that fine-tuning. 



V. THE TOP QUARK EDM 

Now that the top quark has finally been discovered, one can envision some nice experi- 
ments which measure properties of this known particle. Future colliders, such as the NLC, 
can provide many precision measurements of the production cross-section and decay prop- 
erties of the top quark. It is possible that signatures of new physics could arise out of such 
a study. One property of the top quark which has received much attention [Tj| is the pos- 
sibility of measuring its EDM by looking at the decay distributions of the ti pairs. (Other 
CP-violating observables are possible, such as those arising from t — * bW decays, but we 
will make our point only with the top EDM.) It is generally estimated that the top quark 



EDM (d t) can be measured to values as low as O(10 18 )ecm ||19|| . Given the constraints 
which we derived above, we ask if the minimal supersymmetric standard model can yield a 
value for d t this large. 

In the context of supersymmetry, it has been proposed |2(J that a large d t is possible if 



the phase of Af w is of order one. But in Sec. Ill, we showed that Af w is constrained to 
be smaller than or of order the phases which contribute to d n . The EDM of the top is thus 
constrained to be less than a constant times the neutron EDM: 

d t < m detMl 

d n ~^m d detMf 1 j 

where detM| = m^m?, is the determinant of the (down) squark mass-squared matrix, and 
the value of £ depends upon many different SUSY parameters, but is generically of order 1. 
Normalizing d n to the experimental bound, we see that 

detM? d expt 

dt < C A 7 2in % 2 x 10- 21 ecm. (37) 

detM- 10~ 2b e cm 
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In addition to this constraint, we recall that the phase of A t at the EW scale must be less 
than about 1/3, just from the RGE suppression factor r t . Thus, as long as detMj ~ detM^, 
we expect d t to fall about three orders of magnitude below detectability at proposed future 
high energy colliders. 

We can turn this analysis around. If a large top quark EDM is discovered, can it be 
explained in the MSSM? One possibility is that a conspiracy occurs between several large 
phases in the theory to render d n below experimental limits, yet produce a d t detectable at 
high energy colliders. This is equivalent to saying that all the 0(1) coefficients which we 
absorbed into the parameter £ in Eq. (|37|) actually conspire to give £ > 10 3 . As we argued 
in the Introduction, we would not view this as a likely explanation. 

Another possibility to consider is that the top squarks are much lighter than the other 
squarks. For d t to be observable, we would need the determinants in Eq. (|36| ) to have a 
ratio > 10 3 . This is possible, but it too would require some fine-tuning. The large top- 
quark-induced running of the goes in the right direction — the lightest top squark mass 
eigenvalue tends to be smaller than the other quarks. However, i 2 generally tracks fairly 
well with the other squarks, q~L, and thus, we estimate that 

detM] m| 

detM? ~ ^F' (38) 



which means that we would need m~ tl < mgj vlOOO to yield an observable d t . If experiment 
determines that mf > 80 GeV then this condition would imply that the superpartners of the 
light quarks are above 2.5 TeV. This is essentially the heavy squark "solution" to the CP 
violation problem we mentioned in the Introduction, with an additional fine-tuning implied 
by the small ratio Tn^/m^. 

Finally, one could appeal to differences between d t and dd due to effects proportional to 
m 2 /v 2 , which are negligible in dd- To achieve £ of order 10 3 , one again needs a fine-tuned 
conspiracy of couplings. 

Thus we conclude that if a large dt were found, one would probably have to look beyond 
the MSSM for an explanation. 



VI. CONCLUDING REMARKS 

It has long been noted that the phases of soft supersymmetric parameters generically 
lead to an unacceptably large neutron EDM. This fine-tuning problem has slowly become 
less vexing as the theoretical expectations for the squark masses have risen faster than 
the experimental bound on the neutron EDM has fallen. Nevertheless, for squark masses 
below about a TeV, we showed in Sec. V that the phase of B and universal phase of A do 
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not meet the fine-tuning criterion set forth in the Introduction (see Fig 6c). Certainly, if 
super symmetry is discovered at LEP 2 or TeV 33, a fundamental explanation for the absence 
of a neutron EDM would be needed, and any scheme for baryogenesis at the EW scale would 
require that mechanism to leave small effective low energy phases in the soft terms fI5|-|l8|| . 

From the phenomenological point of view, it is tempting to postulate that the soft phases 
are not universal — that the EW phase of A t is large, while the other phases which directly 
contribute to the neutron EDM are small. This would allow interesting signatures of super- 
symmetric CP violation to be visible in top quark physics at future colliders. But we have 
demonstrated by using the renormalization group equations that the imaginary part of A t 
must be less than twice the imaginary part of B, and ^-induced CP-violating observables 
such as the top EDM are thus expected to be unobservably small in almost all minimal 
SUSY models. 

These constraints are particularly important for models of EW baryogenesis which rely 
upon the phase of the stop LR mixing parameter, A t — fitan/3, to generate enough CP 
violation for baryogenesis. Such models must also have sufficiently small \A t — /itan/3| to 



ensure that the phase transition is first order 27 ] . There has also been a recent attempt to 
explain the observed CP violation in the neutral kaon system with zero CKM phase and non- 
zero off-diagonal phases in the general A matrices . If the universal diagonal A parameter 
has a large phase at the GUT scale, it will, as we noted above, give a large contribution to d n 
through a renormalization group induced phase in /i, as well as from a direct contribution. 
One could evade such bounds by insisting that the off-diagonal components of the A matrices 
have a large phase, while the phases of the diagonal A's and of B vanish. Although this 
hypothesis can probably be technically consistent with our fine-tuning criterion (phases 
either zero or large), this scenario strikes us as unnatural. 

Acknowledgements. We would like to thank G. Kane, S. P. Martin, D. Wyler, S. Thomas, 
A. Riotto, and A. Soni for helpful discussions. RG greatly appreciates the hospitality of the 
Brookhaven National Lab HEP Theory Group. 



VII. APPENDIX A 

In this appendix, we provide the details related to our analytic results of Sec. III. It 
is interesting to note that we can use the RGEs for the top Yukawa and gauge coupling 
constants in Eqs. (|9|) and (|TT| ) to write a pseudo-analytic solution to r t = lmAf w /lmAf UT . 
The integral in Eq. (|26| ) can be rewritten as ln(af w /af UT ) — S (c a /6 a ) \n(a^ UT /a^ w ), 
which allows us to write a pseudo-analytic r t in terms of EW and GUT scale quantities (the 
latter of which cannot be found analytically): 
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To place an analytic upper bound on r t , we must place a lower bound on the area 
ItEw T ^ a t(t)dt. We will need the ab = limit of the running of the top Yukawa coupling in 
(I), 

^ = -f(t)a t + 12al (40) 

where f(t) = 2c a a a . While this cannot be solved analytically, we note that as(t) runs down 
with energy and one can show that f(t) will be at its maximum value at the EW scale. 
Thus if we take f(t) to the constant Jew-, we will minimize the running of a t , and Eq ( f40|) 
can be solved analytically to yield the bound 

a,(t) > + - J!) eW.-wA" _ (41) 

\JEW \OC t JEW j ) 

which is valid for 12af w < Jew (larger af w allows the bound on at(t) to reach infinity for 
t < tcur and thus makes the bound useless), which corresponds to tan (3 > 1.3 for m t = 175. 
If we replace a t (t) in the integral above by the RHS of Eq. (0), we can find an analytic 
solution for the lower bound on the area which, for the relevant range of Jew and tEw, can 
be approximated by 



— f'EwtEW — In 




EW S 



EW 



In 



I2a t 



EW 



EW 



(42) 



which yields Eq. (|^) directly. 

For moderate tan/3, we need to include mj effects which mix A t with A^, and a t with 
The coupled differential equations (|17D and (|i~8| ) can be solved analytically only if the 
coefficients, which here are proportional to at and ah, are constants. To obtain bounds on 
the running of A t and Ab, we can break up the range of energy from tsw to tcur into small 
regions where the coefficients are effectively constant, and iteratively evolve from the GUT 
scale down to the weak scale. At each energy tj, the value for A t is given by 

At (t j+1 ) ~ At(t s ) exp (-12^(^)5*) 

( * b{tj) ... [exp {-Ua^t) - exp (-12a t (^)]l , (43) 

Oi t {tj) - 0£ b [tj) J 

provided that T = ab/ (at — &b) is not large. Here 8t = tj — tj+i, which is positive. Iterating 
Eq. p3] ) gives a complicated expression with terms proportional to each of the T(tj)'s. 
However, each of these terms is positive, so that taking T(tj) to its maximum value max- 
imizes the size of the quantity in { }'s in Eq. (|43"D , which is what we need for the case 
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Af UT /Af UT < 0. Once we take T(tj) — > T max , many terms cancel, and we are left with 
(taking 5t — > 0) the upper limit 



Af w < Af UT ex V (- f I2a t {t)dt] 

\ T / r t<3UT \ ( r t<3UT N 

- ) exp - / 12a b (t)dt - exp - / 12a t (t)dt 



6 \ctt — ah. 



(44) 



One can show analytically that T(t) reaches its maximum value at the lowest energy of the 
range, and thus we can replace T max by af w / (af w — a^ w ) . To obtain a simpler bound, 
one can reduce the [ ]'s in Eq. (|44]) to 1 by taking a lower bound on a^(t) to be zero and an 
upper bound on a t (t) to be infinity. Finally one uses the mb — bound on r t obtained in 
Eq. (p8|) for the first term in Eq. ([44]) to obtain the upper bound on r t in Eq. (|30|). 



VIII. APPENDIX B 



In this Appendix we present analytic expressions for the full one-loop SUSY contribution 
to the neutron electric dipole moment, d n . The gluino JTTJ] and chargino |7| contributions 
appear in the literature. While an expression for the neutralino contribution is given by 
Kizukuri and Oshimo |7|, it is written in terms of 4 x 4 complex unitary matrices which 
must be determined numerically. Below we give an expression for this neutralino contri- 
bution solely in terms of the mass matrices (and other MSSM parameters), and a useful 
approximation to that expression, which do not require calculating complex unitary matri- 
ces. 

To find the neutron EDM, we first calculate the EDM of the up and down quarks (d q ) 
from one loop diagrams with photons attached to either (a) an internal boson or (b) an 
internal fermion line. Then the neutron EDM is related to the quark EDM's in the Naive 
Quark Model by d n = {Add — ^«)/3, though recent work has argued that this expression 
overestimates d n if the strange quark carries a large fraction of the neutron and proton spin 
]. The Feynman integrals associated with (a) and (b) are 



[x 



x) 



3 x 
-2 + 2 



In x 

1 — x 



I b (x) 



(1-x) 



1 x a; In a; 

2 + 2 + 1 - x 



(45) 



As we mentioned earlier, all SUSY CP violating effects arise from diagonalizing complex 



mass matrices ]I5]. Gluino loops contribute to the quark EDM d q through the complex 
phase in the left-right mixing elements for up and down squarks: 



-2 m q mglm(A q - iiR, 
d q {9) = ^-Q q ea s j 



'/(I 



ni 



(46) 



go 
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We have averaged over the nearly degenerate squark mass eigenstates for simplicity: m 



and I(x ) = (I(xi) +/(x 2 ))/2. Here Q q e is the charge of quark q, R q = tan/3 (cot (5) 
for q = d (u), and rrig is the gluino mass. (Note that we use Imz for the imaginary part of 
z in this appendix because it is clearer than z in more complicated expressions.) 

The chargino contribution is proportional to the imaginary part of products of elements 
of the matrices U and V which diagonalize the chargino mass matrix. It turns out that one 
can write those products directly in terms of the elements of the chargino mass matrix, so 
that the chargino contribution to d q can be written 



t m q m w lmfi [uoI a ( yi ) + g^/Q] - [ool a (y 2 ) + Q q >I b (y 2 )} 

(47r) 2 \ y 1Lg mi 2/i — 2/2 



^-M«(y») + IT lm{Ar ; ^ m ^' 2 ), 



sin p cos p _ , rrif v z 



r=ct i 



where is the wino mass, to — +1 (—1) for g = d («), j/i = m^/m?/, and I(j/o) = 
(I(yx) + I(y-2))/2. The primed quantities refer to the SU(2) partner, so if q = d, then q' = u 
and r is summed over the set {u,c,t}. Previous expressions for d q (x + ) have neglected the 
squark mixing piece, which is the second term in Eq. (0). This piece is suppressed relative 
to the other contributions by m 2 \ V qr \ 2 /v 2 , which is less than 10~ 4 for q = u or d (but it 
can affect the EDM's of other quarks), though it is interesting that there is a (tiny) direct 
contribution to d n from Inn4 t . 
The neutralino contribution, 

RelAM)*, ^^ _ l m{ A, j R> , ^ alARet ] (4g) 
m qo h,l m qo h,l ) 

arises from the 4x4 complex neutralino mass matrix. The index q = 3 (4) for q = d (it). 
Recall |!T] that the "1" and "2" weak eigenstates are gauginos, and the "3" and "4" weak 



eigenstates are higgsinos which couple to down and up quarks respectively. Thus "34" and 
"43" terms are absent, which will allow us to simplify expressions involving the neutralino 
mass matrix, since that is the position of the complex coefficient //. We have dropped terms 
of order m 2 /v 2 relative to the others. The gauge coefficients an are: 

a q L1 = y/2gtaxiO w (Q q - T q iL ) = V2gttm6 w /Q, (49) 
a£ 2 = V2gT* L , (50) 

and the am are the same as the an with T%l — > T3R = 0. The neutralino phases appear 
through a 4x4 matrix 
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*M = £t/£M«/£0ii, (51) 
1=1 

where [/ diagonalizes the neutralino mass matrix M, and M is the diagonal result. Here 
l\j = I b (xi)5ij, i.e. Ifj is the diagonal matrix of Feynman integrals for the corresponding 
mass eigenvalues in M. In the limit that the I b (xi) are equal, the real part of $^ can simply 
be written 

Re$ w ~ I b (x )ReM hl , (52) 

where I b (x ) = Yh I b (xi)/A. The imaginary part of $ w is more difficult because it vanishes 
in the limit of degenerate neutralino masses (except for the irrelevant "34" and "43" terms) . 
We know that Im/z is the only complex coefficient in the neutralino mass matrix M, so we 
can write 

hn® h i ~ Q M 1mfj,, (53) 

where Qhi is a real matrix to be determined. This allows us to see that d q (x°) is proportional 
to Im/i and lm(A q — fiR q ), just as the other contributions are. It turns out that Im$ w is 
proportional to lm(MM*M) hl , lm(MM*MM*M) hl , and Im(MM* MM* MM* M) hl (except 
for the "34" and "43" pieces). To extract the Im/i dependence, we ignore all terms of higher 
order in Im/x/|/x|, which is a valid approximation for the phases allowed by the experimental 
bound on d n . Then these products (for (h,l) ^ (3,4), (4,3)) simplify as follows: 

lm(MM*M) hl ~ Im^(-l) p (M p R PM 2 R - p ) hl , (54) 

p=0 
4 

lm(MM*MM*M) hl ~ Im/i ^(-lf(M^M^) w , (55) 

p=0 
6 

lm(MM*MM*MM*M) hl ~ Im 1 u^(-l) p (M^PMj|" p ) u , (56) 

p=0 

where M# = ReM and P is a matrix with —1 in the 34 and 43 positions and everywhere 
else (so that Im(//P) = ImM). After some calculation, we obtain an expression for the 
imaginary part of the complex matrix $: 

3 4 ^ s 
lm$ hl ~ Qhilm/j, ~ -lmfi^2(l(x 8 ) - E 1 ^')) x 

1 8=1 j 

E e ijkMfMj £(-lRM£PA4-% - £ ey fc JW?Mj E(-l) P (M^PM^) h/ + 

i,j,fc p=0 p=0 

E^M/E(-1) P (M^PM^)J x (57) 

p=0 
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E etikMfM! (3M S 2 - E M l) - E (SM S 4 - E ^ n 4 ) + 

.i,j,k n j.jM n 

E^M-M/(3M S 6 -E^ 

i,j,k n 

where J2f s means sum over the three members of the set {1, 2, 3, 4} — {s}. Here Mj are the 
mass eigenvalues of M (i.e. the four physical neutralino masses). 

The expression above is completely analytic and exact except for the approximation we 
made in dropping higher order terms in Im/x/|/x|, but it has so many terms that it is not 
that useful. Let us find an approximation to this expression using the information about 
the neutralino mass eigenstates, namely that they are fairly close together and the heaviest 
neutralino is lighter than the squarks (x± <C 1) in most SUSY models. This means that we 
can take a simple linear fit to the Feynman integral by evaluating I b (x) at the lowest and 
highest values of x: 

l\x) ~ K + K x x ~ I b (xi) + S u (x - xi), (58) 

where S41 is the slope 

= />(».)-/><»,) 

X4 — X\ 

and Xj = m^o/rrigg. Thus K x = S41 and K = I b (xi) — 541X1. Note that this approximation 
gives exact values for x\ and X4, and is only off for x 2 and 23 — a rough estimate is that the 
approximation is correct to about 5%. 

If we plug Eq. fl5"8"D into Im$^ in Eq. (0), we see that the K piece vanishes (except 
for the "34" and "43" pieces), and we obtain 



'hi 



lm(MM*M) hl S 41 2 
5 « —2 - -?T lm ^ El -1 ) ( M R PM R )hl- (60) 



m g0 m q0 p=0 



The neutralino contribution to d q is found by plugging Eq. ( [52]) for Re$^ and Eq. (|57|) or 
Eq. (H) for Im$ w into (g|). 

Finally, we want to relate the expressions for the three SUSY contributions to the quark 
EDM in Eqs. (fR]), (0), and ( f48|) in terms of the coefficients k n from Section IV. Using the 
Naive Quark Model, d n = 4/3<i d — l/3<i u , we can write k x = 4:/3k% — l/3k* and 

K = 77^^' < 61 > 

y 1(J zo e cm Lmx 

where x = A u , A^, or /i, and d x q is the contribution to d q from complex quantity x. We 
can see that if we neglect the tiny second term of d q (x + ) m Eq. (0), then k^ u (k^ d ) gets 
contributions only from the gluino and neutralino contributions to d u (da), whereas k£ gets 
contributions from all three of the SUSY contributions to d u and dd- 
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FIG. 2. The ratios of imaginary parts to lmAf UT versus tan/3 with lmAf UT / and 
ImAf UT = 0. The solid line is ImAf w /ImAf UT ; the dashed line is lmB EW /ImAf UT ; and the 
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FIG. 3. Same as Fig. 2 with ImAf UT = ImA? UT . 
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FIG. 5. Scatter plot of (a and d) k^ u , (b and e) k^ d , and (c and f) k£ versus squark mass 
and versus tan (3. Each point represents a solution of the super symmetric parameter space with 
universal scalar and gaugino mass terms which is within other experimental limits. 
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FIG. 6. Scatter plot of d n /10 _25 e-cm versus squark mass and versus tan/3 for (a and d) 
ArgA GUT = ArgA GUT = 0.1 with all other phases zero, (b and e) ArgA GUT = -ArgA GUT = 0.1 
with all other phases zero, and for (c and f) universal phases ArgA GUT = ArgB GUT = 0.1. Each 
point represents a solution of the supersymmetric parameter space with universal scalar and gaug- 
ino mass terms which is within other experimental limits. 
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FIG. 7. Scatter plot of <i Tt /10~ 25 e-cm versus ImAf w /m®. Each point represents a solution of 
the supersymmetric parameter space with universal scalar and gaugino mass terms which is within 
other experimental limits. 
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